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Abstract 

This paper discusses the progress of work to model high-speed supersonic reacting 
flow. The purpose of the work is to improve the state of the art of CFD capabilities 
for predicting the flow in high-speed propulsion systems, particularly combustor flow- 
paths. The program has several components including the development of advanced 
algorithms and models for simulating engine flowpaths as well as a fundamental exper- 
imental and diagnostic development effort to support the formulation and validation 
of the mathematical models. The paper will provide details of current work on exper- 
iments that will provide data for the modeling efforts along with with the associated 
nonintrusive diagnostics used to collect the data from the experimental flowfield. Sim- 
ulation of a recent experiment to partially validate the accuracy of a combustion code 
is also described. 


Introduction 

The design and development of a propulsion system is accomplished with several levels of 
analytic tools, ground-based testing and finally flight. This paper will discuss the progress 
of work to improve the state of the art of CFD capabilities for predicting high-speed re- 
acting flow in engines. The research has several components including the development of 
advanced algorithms and models for simulating engine flowpaths as well as a fundamental 
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experimental and diagnostic development effort to support the formulation and validation 
of the mathematical models. The paper will provide details of current work on experiments 
and the associated non-intrusive diagnostics that will provide data for the modeling efforts. 
This effort will lead to the development of phenomenological models for Reynolds averaged 
Navier-Stokes codes, subgrid scale models for large-eddy simulation techniques, and reduced- 
kinetics models that can then be applied to the design of high-speed propulsion systems. 

Axisymmetric Coaxial Free Jet Experiments and Related Design 

Computations 

Computational fluid dynamics (CFD) methods that employ the Reynolds-averaged Navier- 
Stokes (RANS) equations are widely used in the design and analysis of hypersonic airbreath- 
ing engine flow paths. These methods require models of statistical quantities of the tur- 
bulence in their development which have to be empirically calibrated and validated. In 
particular, new models for turbulent Schmidt and Prandtl number, as well as for turbulence- 
chemistry interactions, are required . 1 While suitable data is available for low-speed flows 
with combustion, sufficient data is still lacking in the area of supersonic combustion. Goyne 
et al. report measurements using particle-imaging velocimetry of mean streamwise velocity 
in a dual-mode scramjet . 2 At the NASA Langley Research Center, several data sets have 
been acquired in a H 2 fueled supersonic combustor using the coherent anti-Stokes Raman 
spectroscopy (CARS) technique 3 and the dual-pump CARS technique . 4,5 The former tech- 
nique gave temperature only while the latter gave both temperature and composition. Data 
included both mean flow and turbulent statistics, although the uncertainty in the latter was 
limited by both instrument precision and number of measurements from which the statistics 
were formed. International work in this area includes measurements in scramjet combustors 
conducted at ONERA (France) and DLR (Germany) using CARS , 6 and other non-intrusive 
techniques. 

Available data arc limited to only a subset of the important variables (which are temper- 
ature, composition, density and velocity) in a limited number of geometries, and turbulence 
data are even rarer. To meet the need for more data, an Interferometric Rayleigh Scattering 
technique for measuring instantaneous velocity is being developed at the NASA Langley 
Research Center to complement the Dual-Pump CARS technique . 7,8 This technique uses 
Rayleigh scattering from one of the CARS laser beams to measure velocity in the same 
instant that CARS measures temperature and composition. Details on the techniques are 
given in a later section. Analysis of streams of such instantaneous measurements allows 
formation of the statistical quantities (means, variances, and covariances) required by mod- 
elers. Experimental configurations to provide suitable flow fields are being developed. An 
axisymetric coaxial free jet was selected since it provides the good optical access required for 
the Rayleigh technique, and symmetry can be taken advantage of to minimize the number of 
required spatial points. In order to form accurate statistics, large numbers of measurements 
are needed at each given point. The use of symmetry allows the total number of measure- 
ments to be made manageable. Another requirement of the flow held is that it should provide 
data that is relevant to both H 2 - and hydrocarbon-fueled hypersonic air-breathing engines. 
The experiment should be capable of providing hows of various Mach numbers to establish 
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Figure 1: The supersonic jet combustor and nozzle, (a) Sectional view (without bolts, tubes, 
etc.), (b) Image during CARS-Rayleigh optical system data acquisition. 


the effects of compressibility. Finally, it is desirable to be able to create both supersonic com- 
busting flows in which the flame is attached to the burner (flame held) and flows in which 
the flame is detached, since both types of flow can occur in hypersonic engine combustors. 
Figure 1 shows a schematic and photo of the coaxial free jet experiment. Figure 1(a) shows 
the burner and nozzle, sectioned along the axis, without bolts, tubes, spark plug and other 
fittings shown; Figure 1(b) is an image of the facility near the nozzle exit during operation, 
showing the supersonic jet of vitiated air and the laser beams of the CARS-Rayleigh system. 
The facility consists of a water-cooled combustion chamber (burner), a silicon carbide nozzle 
(sonic convergent, or supersonic convergent-divergent for M=1.6 and M=2), with an exit 
diameter of 10 mm and a coflow nozzle. Various combinations of gas flows to the burner 
are possible. In one set of experiments, II> or C II 4 fuel, air and 0 2 are reacted to provide 
vitiated air at various temperatures (dependent on flow rates). In these experiments, the 
coflow is of unheated H 2 or CH 4 , and the result (depending on temperature) is a mixing and 
reacting coaxial jet flow. In another set of experiments, the gas flows to the burner are H 2 , 
0 2 or N 2 , and air in such a ratio that the products contain excess H 2 . The coflow is of air 
and the result again is a mixing and reacting coaxial jet flow. Depending on temperature 
and on the co-flow rate, the flame may be held at the annular base region formed between 
the central and coflow nozzle exits where the flow recirculates. Non-reacting cases can also 
be considered in the case where there is no excess H 2 in the center jet and no fuel in the 
co-flow. 
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Figure 2: OH contours showing the sensitivity of the flame location to the presence of angled 
co-flow injection. 

As part of the design effort, a parametric study of the proposed geometry and run con- 
ditions was conducted using CFD to determine the sensitivity of important experimental 
quantities, e.g. temperature, species, density and velocity, to a number of variables includ- 
ing: inner jet temperature, flame holding base size, co-axial jet injection angle, inner jet 
Mach number and sensitivity to computational values of the turbulent Prandtl and Schmidt 
numbers. An example calculation using the Vulcan CFD code 10 is shown in Figure 2, which 
shows the sensitivity of the flame location to the presence of a co-flow. This flow is with a 
Mach 2, hot hydrogen-rich inner jet without and with cold air co-flow. 

Once the experimental geometry was finalized, CFD was used to help define the ex- 
perimental test matrix. This matrix, given in Table 1, includes variations in the inner jet 
Mach number, the temperature of the inner jet, the type of fuel (hydrogen or methane), the 
equivalence ratio of the fuel and air and the location of the fuel and air streams (inner jet or 
co-flow jet). The CFD calculations completed in this step provide a baseline to compare with 
improved models planned for development using data from this and other experiments in 
the program. Figure 3 contains tables and flow images showing the various types of flames 
that are observed in the actual experiment. Two types of images are shown: visible light 
(true color images) with a long exposure time, and false color infrared (IR) images acquired 
in the 8-9 micron (long wave) region at an exposure time of 10 ms. The tables contain 
information pertaining to the state of the flame: no flame, detached flame, flame held at 
the base or at the external coflow boundary, and if the flame holding is marginal, i.e. , at 
the point of extinction. Figure 3(a) pertains to cases with ^-vitiated air in the center jet 
and H 2 in the coflow. The center jet sensible enthalpy is varied by adjusting the flow rates 
from that of Mach 5.5 flight to that of Mach 7, and the exit Mach number of the center jet 
is either 1.0 or 2.0. The coflow is subsonic at the nozzle exit and the equivalence ratio (the 
ratio of H 2 coflow rate to that required to consume all the O 2 in the center jet, 0) is either 


4 


Inner Jet 

Outer Jet 

Mach No. 

Heater Operation 

Mach No. 

Gas 

1,2 

H2+O2+ Air vitiated no 

unreacted H 2 

one T 0 

off 

- 

1,2 

H2+O2+A1T vitiated 

II) rich 

various T a 

off 

- 

1,2 

H2+O2+A1T vitiated 

II) rich 

various T a 

<1 

Air 

1,2 

H2+O2+A1T vitiated 

O2 rich 

various T a 

<1 

ch 4 

1,2 

CH4+O2+A1T vitiated 

O2 rich 

various T a 

<1 

ch 4 

1,2 

H2+O2+A1T vitiated 

O2 rich 

various T a 

<1 

h 2 

1,2 

CH4+O2+A1T vitiated 

O2 rich 

various T 0 

<1 

h 2 


Table 1: The experimental test matrix 
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Figure 3: Flame state for a matrix of supersonic operating conditions at sonic and Mach 
2 pressure matched exit conditions. Mixed infrared (false color) and visible light images 
(two right-most images), (a) Vitiated air center jet with subsonic hydrogen air coflow at an 
overall equivalence ratio (j). ( b) Hydrogen rich center jet with sonic coflow of air. 
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Figure 4: Coaxial jet assembly cross-section 


0.5 or 1. (It is not implied that the coflow H 2 reacts only with the 0 2 in the center jet.) For 
M exit = 1, the flame is held at the exterior boundary in all cases, whereas, for M exit = 2, in 
some cases there is no flame and in others the flame is detached (stands off from the nozzle 
exit), but there are no cases with flame holding. With a detached flame, the trend with 
increasing center jet enthalpy is for the flame to move towards the nozzle exit. Figure 3(b) 
pertains to cases with excess H 2 in the center jet, either 33% or 50% by volume of the jet 
flow being H 2 . The exit static temperature is increased in increments of 100 K, and the exit 
Mach number is either 1.0 or 2.0. The coflow of air is sonic and pressure matched at the 
exit. In some cases, no flame is observed; in others, detached flames or flames attached at 
the base. The trend with increasing T ex n and decreasing M ex n is from no flame to detached 
flame to base-attached flame. Thus, a full range of flame behavior is attained which is ideal 
for model development and code validation. 

Work has just begun to simulate the axisymmetric coaxial free jet experiment utilizing 
both RANS and LES codes. This work has not yet been completed. Work has been com- 
pleted, however, on a coaxial jet mixing experiment that preceded the current work, and this 
experiment has been simulated using RANS codes. Details of the experiment are described 
in two earlier papers. 11,12 Low-density helium, which serves as a simulant of hydrogen fuel, 
was chosen to allow detailed studies of mixing without chemical reaction. Oxygen is added 
to the helium jet as a diagnostic aid for an oxygen flow-tagging technique (RELIEF). Sev- 
eral methods are utilized to characterize the flow field including Schlicren visualization, pitot 
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Figure 5: Schlieren image of coaxial jet mixing (conical extension cap removed) 


pressure, total temperature, and gas sampling probe surveying, and RELIEF velocimetry. A 
schematic of the coaxial jet configuration is shown in Figure 4. The rig consists of a 10 mm 
inner nozzle from which helium, mixed with 5 percent oxygen by volume, is injected at Mach 
1.8 and an outer nozzle 60 mm in diameter from which coflowing air is introduced also at 
Mach 1.8. The velocity ratio between the two jets is 2.25, the convective Mach number is 0.7, 
and the jet exit pressures are matched to one atmosphere. The resulting flow downstream 
of the nozzles can be seen in Figure 5, which shows a Schlieren image of the flowfield. The 
development of the mixing layer between the central helium jet and the air jet can be seen 
along with the shear layer development between the air jet and the surrounding quiescent 
laboratory air into which the air jet exhausts. Shock-expansion wave structure emanating 
outward from the centerbody nozzle lip can also be seen. Inward propagating waves from 
the inner lip, due to the finite thickness of the lip (0.25 mm), can be observed in the air jet 
once they pass through the helium jet. These waves are not visible in the helium jet due to 
the low refractive index within the center jet. A third wave can also be observed emanating 
inward from the outer nozzle lip and traversing both the air jet and the helium jet. The 
detailed data taken from this experiment along with the imaged shock-expansion structure 
provide an excellent case for assessing the ability of a reacting flow code to simulate complex 
mixing processes in a high-speed flow. 
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Figure 6: Experimental setup of the combined CARS - IRS system 


Simultaneous CARS and Interferometric Rayleigh Scattering 

It was recently reported 13 and is summarized here, the combination of a dual-pump coher- 
ent anti-Stokes Raman scattering system with an interferometric Rayleigh scattering system 
(CARS-IRS) to provide time-resolved simultaneous measurement of multiple properties in 
combustion flows. Time-resolved simultaneous measurement of temperature, absolute mole 
fractions of N 2 , 0 2 , and H 2) and two components of velocity in a Hencken burner flame 
were performed to demonstrate the technique. The experimental arrangement of the com- 
bined system is shown in Fig. 6. For the measurement of temperature and the absolute 
mole fractions of N 2 , 0 2 . and H 2 a dual-pump CARS method 4 is used. The system uses 
spectrally-narrow green (injection seeded Nd:YAG at 532 nm) and yellow (NB Dye laser at 
552.9 nm) laser pump beams and one spectrally broad red laser (BB Dye laser at 607 nm) 
beam as the Stokes beam. 

The beams are combined at the focusing point of a spherical lens Lc (focal length of 
410 mm) in a folded BOXCARS geometry to probe Raman transitions of N 2 , 0 2 and H 2 . 
The input beams plus the coherent blue signal beam at 491 nm are collected and collimated 
by another spherical lens with the same focal length as Lc. The three input beams are 
captured in a beam dump while the blue signal beam is passed to a spectrometer. The 
CARS signal, which is a spectrally broad blue beam that contains N 2 , 0 2 and H 2 spectra, is 
analyzed by a spectrometer and recorded by the CCD1 camera. The shape of these spectra 
provides information on the temperature while the relative intensities of these spectra provide 



Figure 7: Simultaneous spectra of CARS (left) and Rayleigh scattering from two viewing 
directions (right). The CARS spectra are both data and fits of theory to the data. The 
Rayleigh spectra are data points connected by lines 


a measure of the relative mole fractions. The spectrum is fit with a theoretical model to 
determine the temperature and mole fractions. 

The velocity measurement is performed simultaneously using an interferometric Rayleigh 
scattering measurement system.' The same pulsed, seeded green laser beam employed in the 
CARS technique is used as a narrow-band light source for the Rayleigh scattering system. 
The receiving optics for Rayleigh scattering are designed to capture the Rayleigh scattered 
light from the green pump beam in the measurement volume while preserving the scattering 
angle information, and to mix it together with the unshifted light of the laser before it 
is passed through the interferometer. In this way more than one component of velocity at 
multiple points can be measured in one interferogram. 7,14 The velocity components measured 
are those that bisect the angles of laser’s incidence and the collection optics. Since two 
directions are collected, two velocity components are measured. The Rayleigh scattered 
light from the two measurement directions are combined with a polarizing beam splitter 
(PBS) and directed through a Fabry-Perot etalon and onto a CCD camera for analysis 
(CCD2). The particular setup used here gives a range of measurable velocities from 0 to ~3 
km/sec up to temperatures of about 2500 K. The CARS spectra and the Rayleigh spectra are 
acquired simultaneously by synchronizing the cameras with the green laser Q-switch pulse 
at 20 Hz. The spectra are subsequently processed, as described in reference 13 to determine 
the temperature, composition and velocity. To demonstrate the method, experiments were 
carried out in a Hencken burner, which provides an adiabatic H 2 -a ir atmospheric pressure 
flame. The flame was stabilized with a co-flow of N 2 . This flame was used because it provided 
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a challenging high-temperature measurement environment while producing a known (near- 
zero) velocity. 

Simultaneous CARS and Rayleigh scattering spectra up to 1610 K are shown in Fig. 7 
(a-c) and in Fig. 7 (d-f), respectively: (a) and (d) are in the co-flow of N 2 , (b) and (e) are in a 
high temperature flow containing N 2 and 0 2l and (c) and (f) are in a high temperature flow 
containing N 2 and a low proportion of 0 2 . The CARS spectral plots show the experimental 
data, the theoretical fits, and also the residual between them. These spectra were used to 
calculate the rotational- vibrational gas temperature and the mole fraction of H 2 , 0 2 , and 
N 2 . 

For the Fabry- Perot spectra shown in Fig. 7, only one fringe order from each spectrum has 
been analyzed for this paper. The narrow shaped peak, more visible in the spectra of higher 
temperature, is the reference laser frequency (no Doppler shift). The broader spectrum is 
the Rayleigh scattered light. Each figure contains two Rayleigh and two reference spectra 
corresponding to the two collecting directions VI (black trace) and V2 (red trace of smaller 
amplitude), differing by 34 degrees. Two-component velocity measurements in the range 
of 300-1600 K showed near-zero velocities (< 30 m/s) and standard deviations of 30-40 
m/s. These errors are about 1% of the dynamic range of this measurement system (3,000 
m/s). Measurement of one velocity component was shown to be possible at up to ~2,400 
K with this system. These accuracies and precisions are within the desired range required 
for the planned supersonic combustion experiments where velocities will be up to 1,500 m/s. 
However, it is anticipated that future improvements in both hardware and software will allow 
these errors to be reduced by a factor of two or more. 

Simulation of Mixing and Reacting Flows 

Improved numerical simulation of high-speed propulsion systems and engine performance 
relies on the development of enhanced codes with an increased capability to model turbulence, 
turbulent mixing, and kinetics. A significant amount of the design effort utilizes Reynolds 
averaged Navier-Stokes (RANS) codes for flowpath prediction. RANS design codes utilize 
phenomenological models to describe the turbulence field, fuel-air mixing, kinetics, and the 
interaction of turbulence and kinetics. There is an increased reliance on the use of large 
eddy simulation (LES) for the simulation of high-speed reacting flows. 16 LES methods 
calculate the large scale features of the flowfield and require modeling only to describe 
the small scale features of the flow utilizing a subgrid scale model. Subgrid scale models 
using filtered density functions (FDF) have shown great promise in improving the state-of- 
the-art in combustor modeling. 1 ' The FDF method is the counterpart of the probability 
density function (PDF) methodology, which has been long-established in RANS calculations 
of turbulent combustion. 18 The axisymmetric version of the SPARK RANS code was used 
to simulate the flowfield in the helium/oxygen center nozzle and the outer air nozzle of the 
coaxial jet mixing experiment described in the previous experiment. Details of the code are 
given in reference 19. The analysis of the experiment was begun by first solving for the 
flowfield in the center and outer nozzles. Results obtained at the end of each nozzle were 
then used to specify the supersonic inflow conditions for the downstream domain beyond the 
nozzles where mixing of the jets occurred. Details regarding the nozzle calculations and the 
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Figure 10: Comparison of RELIEF 



Figure 9: Comparison of pitot pressure data 
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mixing region downstream of the nozzles are given in reference 20. 

A comparison of the measured helium mass fraction data with the simulation results at 
several stations downstream of the nozzles is given in Figure 8. Agreement between the 
simulation and the data is very good at each station. The code somewhat overpredicts the 
mixing near the centerline at the x = 0.12 m station, although the prediction improves with 
increasing radial distance. A comparison of measured pitot pressures with the simulation 
is shown in Figure 9. Agreement is good in the region of the air coflowing jet, but the 
simulation somewhat overpredicts the pitot pressure in the helium-air mixing region. The 
comparison with the experimental data differs at large radial distances greater than 0.025 m 
as the code does not consider the effects of the laboratory air entrained by the coaxial air jet. 
The RELIEF streamwise velocity data is compared with the simulation in Figure 10. The 
prediction agrees well with the data at the first three stations and slightly overpredicts the 
data at the remaining stations near the centerline. The simulation somewhat underpredicts 
the the velocity at the final three stations in the mixing region between the helium and air 
coflowing jets in agreement with the pitot pressure results. 

Concluding Remarks 

This paper has described work to develop phenomenological models for Reynolds aver- 
aged Navier-Stokes codes, subgrid scale models for use in large-eddy simulation and reduced- 
kinetics mechanisms to model hydrogen-air and hydrocarbon-air reaction in propulsion ap- 
plications. Fundamental experiments are being performed to provide data that will be used 
in the development and refinement of these models. Experimental data is extracted from the 
experiments using nonintrusive diagnostics that allow accurate simultaneous measurement 
of temperature, species, and up to three components of velocity in supersonic flow without 
changing the character of the flowheld. Once the databases are in hand, the data will be 
analyzed using a response surface methodology 21-23 that provides an efficient means of de- 
termining critical parameters in a chosen model. The models will then be incorporated into 
combustion codes used in engine flowpath analysis and design. 
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